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Abstract. Exploiting symmetry in Grobner basis computations is difficult 
when the symmetry takes the form of a group acting by automorphisms on 
monomials in finitely many variables. This is largely due to the fact that the 
group elements, being invertible, cannot preserve a term order. By contrast, 
inspired by work of Aschenbrenner and Hillar, we introduce the concept of 
equivariant Grobner basis in a setting where a monoid acts by homomorphisms 
on monomials in potentially infinitely many variables. We require that the 
action be compatible with a term order, and under some further assumptions 
derive a Buchberger-type algorithm for computing equivariant Grobner bases. 

Using this algorithm and the monoid of strictly increasing functions N — > N 
we prove that the kernel of the ring homomorphism 

I i, 3 6 N,i > j] -> K[si,tj | i £ N], y i3 h-> SjSj + Ut 3 

is generated by two types of polynomials: off-diagonal 3 X 3-minors and pen- 
tads. This confirms a conjecture by Drton, Sturmfels, and Sullivant on the 
Gaussian two-factor model from algebraic statistics. 



1. Introduction and results 

Equivariant Grobner bases. Algebraic varieties arising from applications often 
have many symmetries. When analysing such varieties with tools from compu- 
tational algebra, it is desirable to do so in an equivariant manner, that is, while 
keeping track of those symmetries, and if possible exploiting them. The notion of 
Grobner basis, which lies at the heart of computational algebra, depends heavily on 
choices of coordinates and of a term order, which is a well-order on monomials in the 
coordinates. It is therefore natural, at least from a computational point of view, to 
study symmetries of ideals that preserve both the coordinates and the term order. 
Now the term symmetry is usually reserved for certain invertible maps, but it is 
easy to see that an invertible map cannot preserve a well-order; see Remark |2.1| 
Hence we are led to relax the condition that symmetries be invertible. On the other 
hand, if a non- invertible map is to preserve the restriction of the term order to the 
set of coordinates, then that set better be infinite, in contrast with the usual set-up 
in computational commutative algebra. 

In fact, there is another, more compelling reason for allowing infinitely many vari- 
ables: many varieties from applications come in infinite families, and it is convenient 
to pass to a suitable limit. For example, the variety of symmetric n x n-matrices 
of rank 2 has a well-defined projective limit for n tending to infinity, and so does 
the closely related two-factor model that we study in this paper. In both cases, the 
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limit is not only stable under the union of all symmetric groups S n simultaneously 
permuting rows and columns, but also under the monoid Inc(N) of all strictly in- 
creasing maps from N to itself. And while the union of the symmetric groups does 
not preserve any term order, the monoid Inc(N) does preserve such an order; this 
fundamental observation allows us to do computations in Section [3] 

This discussion leads to the following set-up, which we believe will have ap- 
plications to numerous other problems. Let X be a potentially infinite set, whose 
elements we call variables. The free commutative monoid generated by X is denoted 
Mon; its elements are called monomials. Suppose that we have 

EGB1. a term order, i.e., a well-order < on Mon such that m < m' => mm" < 
m'm" for all m, m' , m" € Mon; and 

EGB2. a monoid G, i.e., a (typically non-commutative) semigroup with identity, 
acting on Mon by means of monoid homomorphisms Mon — > Mon preserv- 
ing the strict order: nl = 1, 7r(mrn') = (7rm)(7rra'), and m < m' =>■ irm < 
■Km' for all 7r e G, m, m! G Mon. 

Example 1.1. The setting that Aschenbrenner and Hillar study in \T\ fits into this 
framework, and indeed inspired our set-up. There X = {xi,X2, ■ ■ ■} and G is the 
monoid Inc(N) of all increasing maps 7r : N — > N acting on X by n%i = x n n-\ and on 
Mon by multiplicativity. As a term order one can choose the lexicographic order 
with Xi > Xj if i > j. Aschenbrenner and Hillar have turned their proof of finite 
generation of Inc(N)-stable ideals in K[xi,X2, ■ ■ ■] into an algorithm; see [2]. 

Let K be a field and let K[X] — KMon be the polynomial iT-algebra in the 
variables X, or, equivalently, the monoid if-algebra of Mon. Then G acts naturally 
on K [X] by means of homomorphisms. A G-orbit is a set of the form Gz = {ttz | 
7r € G}, where z is in a set on which G acts. The fact that G acts by monoid 
homomorphisms on Mon implies that the ideal generated by a union of G-orbits in 
K[X] is automatically G-stable, that is, closed under the action of G. 

We use the notation lm(/) for the leading monomial of /, i.e., the <-largest 
monomial having non-zero coefficient in /. By the requirement that G preserve the 
order, we have lm(7r/) = it lm(/). Given an ideal / of -f^[AT], lm(J) is an ideal in the 
monoid Mon, that is, lm(/) is closed under multiplication with any element from 
Mon. If I is G-stable, then so is lm(7). 

Definition 1.2 (Equivariant Grobner basis). Let / be a G-stable ideal in i^[A]. A 
G-Grdbner basis of / is a subset B oi I for which lm(G_B)(= {hn(7r&) | b £ B, ir g 
G}) generates the ideal \m(I) in Mon. If G is fixed in the context, then we also call 
B an equivariant Grobner basis. If G = {1}, then we call B an ordinary Grobner 
basis. 

Remark 1.3. At MEGA 2009, Viktor Levandovskyy pointed out to the second 
author that our equivariant Grobner bases are a special case of Grobner S'-bases 
in the sense of [B], which were invented for analysing certain two-sided ideals in 
free associative algebras. The focus of the present article is on getting exactly the 
right set-up for doing machine computations of equivariant Grobner bases in the 
commutative setting. 



It is easy to see that if B is a G-Grobner basis of /, then GB generates / as an 
ideal; see Lemma 2.2 
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Example 1.4. Let X = {y^ | i,j 6 N}, let k be a natural number, and let / be 
the ideal of all polynomials in the that vanish on all N x N-matrices y of rank 
at most k. Order the variables yij lexicographically by the pair where i is the 

most significant index; so for instance > j/2.6 > J/2,4 > 2/1, io- The corresponding 
lexicographic order on monomials in the t/y is a well-order. Let G :— Inc(N) xlnc(N) 
act on X by (jr, o~)yij = y^(i)^(j)\ this action preserves the strict order. The G-orbit 
of the determinant D of the matrix (yij)i,j=i,., ,,k+i consists of all (k + 1) x (fc + 1)- 
minors of y, which by the results of [TT] form a Grobner basis of the ideal /. As a 
consequence, {D} is a G-Grobner basis of I. 

A G-stable ideal need not have a finite G-Grobner basis. Indeed, if one requires 
that every G-stable ideal / in K [X] has a finite G-Grobner basis, then this must in 
particular be true for monomial ideals. This implies that Mon does not have infinite 
antichains relative to the partial order on Mon defined by m -< m! 3ir 6 G : 
7Tm\m' . Observe that this is, indeed, a partial order: transitivity is straightforward, 
and antisymmetry follows from the fact that itm\m irm < m, while on the other 



hand nm > m for all 7r, m; see Remark 2.1 Conversely, if (Mon, X) does not have 



infinite anti-chains, then every ideal has a finite G-Grobner basis. This is the case 
in the set-up of Example |1.1| which is generalised in [10] : there equi variant Grobner 
bases are called monoidal Grobner bases. 

Remark 1.5. We have not yet really used that Mon is the free commutative 
monoid generated by X. So far, we could have taken Mon any commutative monoid 
equipped with EGB1 and EGB2. This viewpoint, and a generalisation thereof, is 
adopted in |10j . However, for doing computations we need that Mon has more 
structure; see conditions EGB3 and EGB4 below. This is why we have restricted 
ourselves to free monoids Mon. 



In the polynomial ring of Example[L4|the set {2/122/21, 2/122/232/31, 2/122/232/342/41, ■ ■ •} 
is an infinite ^-antichain of monomials, hence the Inc(N)-stable ideal generated by 
it does not have a finite Inc(N)-Gr6bner basis. But even in such a setting where 
not all G-stable ideals have finite G-Grobner bases, ideals of interesting G-stable 
varieties may still have such bases. We will derive an algorithm for computing 
equivariant Grobner bases under the following two additional assumptions: 

EGB3. for all ir € G and m, m! £ Mon we have lcm(7rm, irm') — -k lcm(m, m'); and 
EGB4. for all f,h € K[X] the set Gf x Gh is the union of a finite number of G- 

orbits (where G acts diagonally on K[X] x K[X}), and generators of these 

orbits can be computed effectively. 
Note that EGB3 is automatically satisfied if G stabilises the set X of variables. 
Although this is the only setting that we will need for the application to the two- 
factor model, future applications may need the greater generality where X is not 
G-stable. Condition EGB4 is of computational importance, as will become clear in 
Section|2] There we also show in Examples 2.6 and 2.7 that these requirements are 
not redundant. 

Theorem 1.6. Under conditions EGB1, EGB2, EGB3, and EGB4 there exists an 
algorithm that takes a finite subset B of K[X] as input and that returns a finite 
G-Grobner basis of the ideal generated by B, provided that it terminates. 

In Scction[2]we derive this algorithm, and in Section[3]we apply it to a conjecture 
concerning a statistical model to be discussed now. 
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The Gaussian two-factor model. The Gaussian k-factor model with n observed 
variables consists of all covariance matrices of n jointly Gaussian random variables 
Xi, . . . ,X n , the observed variables, consistent with the hypothesis that there exist 
k further variables Z\, . . . , Z k , the hidden variables, such that the joint distribution 
of the Xi and the Zj is Gaussian and such that the Xi are pairwise independent 
given all Zj. This set of covariance matrices turns out to be 

F k , n ■■= {D + SS T | D G M„(R) diagonal and positive definite, and S G M nife (R)}, 

where M n ^{M) is the space of real n x /c-matrices, and M n (M) is the space of real 
n x n-matrices. In 7] this model is studied from an algebraic point of view. In 
particular, the ideal of polynomials vanishing on F k , n is determined for k = 2, 3 
and n < 9. The case where k = 1 had already been done in [4]. The authors of [7] 
pose some very intriguing finiteness questions. In particular, one might hope that 
for fixed k the ideal of Fk n stabilises, as n grows, modulo its natural symmetries 
coming from simultaneously permuting rows and columns. For k = 1 this is indeed 
the case, and for arbitrary k it is true in a weaker, set-theoretic sense [5]. In this 
paper we prove that the ideals of i<2,n stabilise at n = 6. To state our theorem we 
denote by ytj the coordinates on the space of symmetric n x n-matrices; we will 
identify y^ with j/y. Recall from [7] that the ideal of F 2j5 is generated by a single 
polynomial 

P '■= Yq S S n ( 7r )^(l)^(2)2/7r(2) : 7r(3)2/7r(3) : 7r(4)?/ 7 r(4),7r(5)y7r(5),7r(l), 

irSSym(5) 

called the pentad. The normalisation factor is important only because it ensures 
that all coefficients are ±1 — indeed, the stabiliser in Sym(5) of each monomial in the 
pentad is the dihedral group of order 10. We consider P an element of Z[yy | i > j]. 
The ideal of F2.6 contains another type of equation: the off-diagonal minor 

M := det(y[{4, 5, 6}, {1, 2, 3}]) G Zfefc | i > j] 

the determinant of the square submatrix of y sitting in the lower left corner of y. 
If / is any polynomial in Rfyy | i > j] that vanishes on i^.n and if we regard / 
as an element of Rkyij | i > j][yn, . . . ,y nn ], then each of the coefficients of the 
monomials in the diagonal variables yu is a polynomial in the off-diagonal variables 
that vanishes on i*2 n , as well. Therefore the following theorem settles the conjecture 
of Drton, Sturmfels, and Sullivant, that pentads and off-diagonal minors generate 
the ideal of F 2 ^„ for all n; see [7J Conjecture 26]. 

Theorem 1.7. For any field K and any natural number n > 6 the kernel I n (K) 
of the homomorphism Ktyij \ 1 < j < i < n] — > K[s\, . . . , s n , t\, . . . , t n ] determined 
by Vij !- > SiSj + Utj is generated, as an ideal, by the orbits of P and M under the 
symmetric group Sym(n). 

Remark 1.8. In [S] it is proved that F 2y7l equals the set of all positive definite 



matrices with the property that every principal 6 x 6-minor lies in -^2,6- Theorem 1.7 
implies an analogous statement for the Zariski closures of -F2,n and F2 6 ■ 



We sketch the proof of Theorem 1.7 which appears in Section [3j We put a 



suitable elimination order on the monomials in y^, i,j G N, i > j, and report on 
a computation that yields a finite Inc(N)-Gr6bner basis for the determinantal ideal 
generated by all 3 x 3-minors of y. Intersecting this Inc(N)-Gr6bner basis with the 
ring in the off-diagonal matrix entries gives Theorem |1.7| 
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2. AN ALGORITHM FOR EQUIVARIANT GROBNER BASES 

We retain the setting of the introduction: X is a potentially infinite set and 
Mon is the free commutative monoid generated by X, equipped with a term order 
(EGB1) preserved by the action of a monoid G (EGB2) which also preserves least 
common multiples (EGB3). Condition EGB4 will be needed only later. 

Remark 2.1. Note that G acts by injective maps on Mon by EGB2. It is essential 
that we allow 67 to be a monoid rather than a group. Indeed, the image of G in the 
monoid of injective maps Mon — > Mon contains no other invertible elements than 
the identity: If it G G then irm > m since otherwise m > nm > ir 2 m > . . . would 
be an infinite strictly decreasing chain. But then if (the image of) 7r is invertible, 
we have irm > m > ir~ 1 m > ir~ 2 m > . . ., another infinite decreasing chain. 

We set out to translate familiar notions from the setting of ordinary Grobner 
bases to our equivariant setting. In what follows the coeffient in / of lm(/), the 
leading coefficient, is denoted lc(/), and lt(/) = lc(/)lm(/) is the leading term of 
/• 

Lemma 2.2. /// is G-stable and B is a G-Grdbner basis of I, then GB = {nb | 
tt G G, b G B} generates the ideal I . 

Proof. If not, then take an / G / \ (GB) with lm(/) minimal. Take b G B and 
tt G G with lm(7r6)|lm(/). Subtracting (lt(/)/lt(7r6))7rb from / yields an element in 
I \ (GB) with leading term strictly smaller than that of /, a contradiction. □ 

Algorithm 2.3 (Equivariant remainder). Given / G K[X] and B C K [X], proceed 
as follows: if 7rlm(6)|lm(/) for some tt G G and b G B, then subtract the multiple 
(It (/) /It (?r&) ) 7r6 of irb from /, so as to lower the latter's leading monomial. Do 
this until no such pair (7r, 6) exists anymore. The resulting polynomial is called a 
G-remainder (or an equivariant remainder, if G is fixed) of / modulo B. 

This procedure is non-deterministic, but necessarily finishes after a finite number 
of steps, since < is a well-order. Any potential outcome is called an equivariant 
remainder of / modulo B. 

Definition 2.4 (Equivariant S-polynomials) . Consider two polynomials b ,bi with 
leading monomials m 0l mi, respectively. Let H be a, set of pairs (cjo,0i) G G x G 
for which G6 x Gb\ = |J, q {(ttoo&o, ttcti&i) | tt € G}. For every element 
((To,cri) G H we consider the ordinary S-polynomial 

1cm (cr mo, eri mi) lcm(cr m , dmi) 

S(er Oo,erit>i) := lc(6i) cr b - lc(6 ) crih- 

tqitiq o\m\ 

The set {S(o~obo, cti&i) | (cto,<ti) G H} is called a complete set of equivariant S- 
polynomials for bo,bi. It depends on the choice of H. Under condition EGB4, H 
can be chosen finite. 
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Theorem 2.5 (Equivariant Buchbcrger criterion). Under the assumptions EGB1, 
EGB2, and EGB3, let B be a subset of K[X] such that for all 6 , b\ G B there exists 
a complete set of S-polynomials each of which has as a G-remainder modulo B. 
Then B is a G-Grdbner basis of the ideal generated by GB. 

We will first prove this for ordinary Grobner bases, and from that deduce the 
theorem for equivariant Grobner bases. The proof for the ordinary case is identical 
to the proof in the case of finitely many variables. We include it for completeness, 
and also because we have no reference where the result is stated for infinitely many 
variables. 



Proof of Theorem 2.5 in case G = {1}. We may and will assume that all elements 
of B are monic. Let / denote the ideal generated by B. If lm(£?) does not generate 
the ideal lm(7) in Mon then there exists a polynomial of the form 

beB 

with only finitely many of the fb non-zero, for which lm(/) is not in the ideal 
generated by lm(_B). We may choose the expression above such that first, the 
maximum m of lm(/^6) over all b for which f b is non-zero is minimal and second, the 
number of b with lm(/b6) = m is also minimal. The maximum is then attained for 
at least two values bo, b\ of &, because otherwise m would be the leading monomial 
of /. Write m; := lm(fei) for z = 0,1, and let t$,ti be such that lcm(mo,mi) = 
torno = t\m\ . Now m = lm(/h )mo = lm(/f, 1 )mi is a multiple of both mo and mi, 
and therefore \m(f ba ) is divisible by t ; set 

lt(Ao) 



.4 



to 

Next consider 

S := S{b ,bi) = t b - tib\, 
where we have used that bo and b\ are monic. As is a remainder of S modulo B 
by assumption, we can write S as a sum X^bes s bb with only finitely many non-zero 
terms that moreover satisfy lm(sf,6) < lm(S') < lcm(mo, mi) for all b. Then we may 
rewrite / as 

/ = / - A(S - £ s b b) = + fb + fb)b 

b b 

where = Ast, and 



-lt(/ 6o ) if 6 = 6 , 

//; = { lc(/ 6o )lm(/ fel ) if 6 = 6i, 
otherwise. 



For all b e B we have 



lm((As b )b) = lm(As b b) < lm fr) l C m(m ,mi) 

to 

= lm(/6 ) TO o = m, 

so for all b we have lm(/^6) < m. Moreover, lm((/f, + f b )bo) is strictly smaller 
than m. Finally, lm(/^ bi) = m. We conclude that either max & lm((/ b + f' b + f b )b) 
is strictly smaller than m, or else the number of b for which it equals m is smaller 
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than the number of 6 for which lm(/(,6) equals m. This contradicts the minimality 
of the expression chosen above. □ 



Proof of Theorem \2.5\ using the ordinary Buchberger criterion. We prove that GB 
is an ordinary Grobner basis of the ideal that it generates. By the ordinary Buch- 
berger criterion it suffices to verify that for all 60,61 £ B and 7To,7Ti € G the 
S-polynomial 5(77060, 7Ti&i) has as a remainder modulo GB. By assumption there 
exists a triple (co, 01, 772) for which (77060, tti&i) = (7720060, 772CT161) and for which 
<S(<To&o, 0161) has as a G-remainder modulo B, which means that it has as an or- 
dinary remainder modulo GB. Since G preserves least common multiples (EGB3), 
we have 

5(772(70^0, 7T201&l) = ^2S(a Q b , Ol&i), 

and applying 77 2 to the entire reduction of S(<Tob , CT161) to modulo GB yields a 
reduction of 5(77060,77161) to 0, as claimed. □ 

The following example shows that EGB3 is not a redundant assumption in The- 
orem [53] 

Example 2.6. Suppose that X = {x, y, z\, Z2, ■ ■ ■} and that the monoid G is gen- 
erated by Inc(N) acting by 772^ = z v i and trivially on x, y, together with a single ho- 
momorphism a : Mon — > Mon determined by ax = x, ay — XZi, and azi = Zi + \ for 
all i. Then G preserves the lexicographic order on Mon for which z^ + i > Zi > y > x 
for all i. Now consider the set B = {y + 1}. We have 

G(y + 1) x G(y + 1) = G(y + 1, y + 1) U G(y + 1, xz 1 + 1) U G(xz 1 + 1, y + 1), 



so we may take H from Definition 2.4 equal to {(1, 1), (1, a), (cr, 1)}. The S-polyno- 
mial 5(y+l, 7/+1) is zero, and the S-polynomials S(y+1, xzi + 1) and S(xz 1 +1, y+l) 
reduce to zero modulo y + l and a(y + l) = xzi + l. Hence we have a complete set of 
S-polynomials oi y + l with itself that all G-reduce to zero modulo B. Nevertheless, 
B is not a G-Grobner basis of the G-stable ideal that it generates, since that ideal 
also contains S(xz\ + 1, XZ2 + 1) = Zi — Zi, which does not G-reduce to zero modulo 
B. 

Here is an example where EGB4 is not fulfilled. 

Example 2.7. Let X = {x, y}, let G be the multiplicative monoid of the positive 
integers, where m acts by x <— > x m and y <— > y m . Now 

GxxGy = {{x\y 3 ) \ i,j e Z >0 }, 

while the diagonal G-orbit of the pair (x l ,y J ) equals {(x al , j/ aj ) | a € Z>o}. Hence 
Gx x Gy is not the union of finitely many G-orbits. 

However, if we do assume EGB4, then every pair (60,61) has a finite and com- 
putable complete set of ^-polynomials and have the following theoretical algorithm, 
alluded to in Theorem [O] We do not claim that it terminates, but if it does, then 



it returns a finite equi variant Grobner basis by Theorem 2.5 

Algorithm 2.8 (Equivariant Buchberger algorithm). 
Input: a finite subset B of ^[^]- 

Output (assuming termination): a finite equivariant Grobner basis of the 

ideal generated by GB. 
Procedure: 
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(1) P := B x B; 

(2) while P ± do 

(a) choose (&o, &i) & P and set P ;= P \ {(bo, 6i)}; 

(b) compute a finite complete set S of equivariant S-polynomials for 

(c) while S ^ do 

(i) choose / € S and set S := S \ {/}; 

(ii) compute a G-remainder r of / modulo -B; 

(iii) if r ^ then set B := B U {r} and P := PU(B x r); 

(3) return _B. 

Note the order in which £? and P are updated: one needs to add (r, r) to P, as 
well. The proof of correctness of this algorithm is straightforward and omitted. 

Remark 2.9. If the partial order ^ on Mon defined in the introduction does not 
admit infinite antichains, then the equivariant Buchberger algorithm always termi- 
nates. Indeed, suppose that the algorithm would not terminate, and let r±, r<i, r^, . . . 
be the sequence of remainders added consecutively to B. Then for all i < j we have 
lm(ri) ^ lm(fj) since rj is G-reduced modulo rj. Using the fact that decreasing 
^-chains are finite, one finds an infinite subsequence r,-j , r^ 2 , . . . with i\ < i% < . . . 
such that lm(r ia ) ^ lm(r ib ) holds not only for a < b but also for a > b. This 
sequence contradicts the assumption that ^ does not have infinite antichains. 

3. AN EQUIVARIANT GROBNER BASIS FOR THE TWO-FACTOR MODEL 



Theorem 1.7 will follow from the following result. Let X = {yij | i,jgN,!> j} 
be a set of variables representing the entries of a symmetric matrix. We consider 
the lexicographic monomial order on Mon in which the diagonal variables yu are 
larger than all variables with i > j, and apart from that > yi* y if and only 
if i > i' or i = i' and j > j' . So for instance we have 

V2,2 > J/1,1 > J/5,2 > 2/4,3- 

Note that this monomial order is compatible with the action of the monoid Inc(N) 
of all increasing maps N —> N. For any polynomial p £ K[X] let l(p) denote the 
largest index of p, i.e., the largest index appearing in any of the variables in any of 
the monomials of p. 

Theorem 3.1. For any field K , let Jjq(K) be the ideal in K[X] generated by all 
3 x 3-minors of the matrix y (recall that we identify y^ for j < i with yij). Relative 
to the monomial order < the ideal Jjq(K) has an Inc (N)-Grobner basis B consisting 
of 42 polynomials. The intersection B D K[yij | i > j] is an Inc(N) -Grobner basis 
of Jfq(K) n K[yij | % > j] consisting of 20 polynomials. The largest indices and the 
degrees of the elements in these bases are summarised in Table pO} 



Remark 3.2. The polynomial with largest index 5 in the Inc(N)-Gr6bner basis 
Br\K[yij | i > j] is the pentad P. The five degree-3 polynomials with largest index 
6 in that Grobner basis form the Sym(N)-orbit of the off-diagonal minor M. All 
19 remaining polynomials are already in the Inc(N)-stable ideal generated by these 
polynomials; this latter statement also follows from the result in [7] that at least up 
to n = 9 the ideal of the two-factor model is generated by pentads and off-diagonal 
minors. The complete Inc(N)-Gr6bner bases of J^(K) can be downloaded from the 
second author's website. 
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Kp) 


3 4 5 6 7 8 9 


#P£B 
degrees 

# P eBn K[ yij \i>j] 

degrees 


1 6 11 10 8 5 1 
3 1 3 6 S 10 ^ 1 3 5 5 5 5 s 5 5 5 1 
1 5 8 5 1 
5 1 3 5 5 s 5 5 5 1 



Table 1. Degrees of polynomials in the Inc(N)-Gr6bner basis B of 
Jfi(K), grouped according to largest index. The first row records 
the largest index, the second row the number of polynomials in B 
with that largest index, the third row their degrees with multiplic- 
ities written as exponents, the fourth row counts the number of 
polynomials containing only variables y^ with i > j, and the fifth 
row records their degrees. 



Remark 3.3. A Grobner basis of the ideal of the two-factor model F 2 , n relative to 
circular term orders was already found in |13j . The proof involves general techniques 
for determining the ideal of secant varieties, especially of toric varieties; see also [T^] . 
The Grobner basis found there, however, does not stabilise as n grows — and indeed, 
circular term orders are not compatible with the action of Inc(N). It would be 
interesting to find a direct translation between Sullivant's Grobner basis and ours. 



Theorem |3. 1| implies Theorem 1.7 



Proof of Theorem \1.7\ It is well known that the (k + 1) X (k + l)-minors of the 
symmetric matrix {llij)ij=i,...,n generate the ideal of all polynomials vanishing on 
all rank-fc matrices (for a recent combinatorial proof of this fact, see [321 Example 
4.12]; in characteristic this fact is known as the Second Fundamental Theorem for 



the orthogonal group). Hence the ideal I n {K) of Theorem 1.7 is the intersection 



of the ideal J n generated by the 3 x 3-minors of (yij)i l j=x t ... t n with the ring K[y v 



i > j]. Theorem 3.1 implies that one obtains a Grobner basis of J„, relative to the 
restriction of the monomial order on K[yij \ i,j £ N, i > j] to K[yij \ l<j<i<n] 
by applying all increasing maps {1, . . . , l(p)} —> {1, . . . , n} to all p € B n K\y^ \ 
i > j] with l(p) < n. Such an increasing map can be extended to an element of 



Sym(n), and Remark 3.2 concludes the proof. □ 



We conclude with some remarks on the computation that proved Theorem |3.1| 
First we need to verify EGB4. 

Lemma 3.4. For all b a ,bi e K[y i:j \ z, j e N,i > j] the set (Inc(N)6 ) x (Inc(N)&i) 
is the union of a finite number of Inc(N)- orbits. 

Proof. Consider all pairs (So, Si) of sets So, Si C N with \Si\ = 1(h) for which 
Sq U Si is an interval of the form {1, . . . , k} for some k, which is then at most 
l(bo) + l(bi). Note that there are only finitely many such pairs (So, Si). For each 
such pair let (ttq, tti) be a pair of elements of Inc(N) such that ni maps {1, . . . , 1(h)} 
onto Si', it is irrelevant how n acts on the rest of N. Then we have 

Inc(N)6 X Inc(N)6i = [j Inc(N)(7r 6 , nih), 

(So, Si) 

where the union is over all pairs (So, Si) as above. □ 
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Computational proof of Theorem \3.1\ The 42 polynomials of B were constructed 
by computing a Grobner basis for Jg(Q) with Singular |9j and retaining only 
those polynomials p for which the set of indices occurring in their variables form an 
interval of the form {1, . . . , k} with k < 9. All elements of B are monic and have 
integral coefficients (in fact, equal to ±1 except for the 3 x 3-minor with largest 
index 3, which has a coefficient 2). By the equivariant Buchberger criterion and 
the proof of Lemma [3. 4| we need only Inc(N)-reduce modulo B all S-polynomials 
of pairs (7To&o, ""i^i) with bo,bi £ B and ^ : {1, . . . , l(bi)} — > N increasing and 
such that im7To U im7Ti = {1, . . . , k} for some k. For instance, for bo = b\ = b 
equal to the polynomial in B with largest index 9, we having to Inc(N)-reduce 
5(7To6, 7Ti6) modulo B for all increasing maps 7ro,7ri : {1,...,9} — » {1, ...,18} 
whose image union is an interval {1, . . . , k}. However, if k — 17 or k = 18, then 
7To& and 7Ti6 turn out to have leading monomials with gcd 1, so these cases can be 
skipped. This reduces the theorem to a finite computation involving polynomials 
with largest indices up to 16, which we have implemented directly in C. Finally, to 
deduce the result for all base fields — and to speed up the computation — we used 
the following argument. Since Inc(N)-B n K [yij | 1 < j < i < n] is a subset of the 
ideal of 3 x 3-minors, it is a Grobner basis if and only if the ideal generated by 
\m(B) has the same Hilbert series as the ideal generated by 3 x 3-minors. Since 
this Hilbert series is known and does not depend on the field [3], we may do all 
our computations over one field and conclude that it holds over all fields. We have 
verified the equivariant Buchberger criterion over F2, which made the computation 
slightly faster than working over Q. □ 
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